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ABSTRACT 



Aims. Within the central parsec of the Galaxy, several tens of young stars orbiting a central supermassive black hole are observed. 
A subset of these stars forms a coherently rotating disc. Other observations reveal a massive molecular torus which lies at a radius 
~ 1.5pc from the centre. In this paper we consider the gravitational influence of the molecular torus upon the stars of the stellar disc. 
Methods. We derive an analytical formula for the rate of precession of individual stellar orbits and we show that it is highly sensitive 
upon the orbital semi-major axis and inclination with respect to the plane of the torus as well as on the mass of the torus. 
Results. Assuming that both the stellar disc and the molecular torus are stable on the time-scale > 6Myr, we constrain the mass of 
the torus and its inclination with respect to the young stellar disc. We further suggest that all young stars observed in the Galactic 
Centre may have a common origin in a single coherently rotating structure with an opening angle < 5°, which was partially destroyed 
(warped) during its lifetime by the gravitational influence of the molecular torus. 

Key words, stellar dynamics — Galaxy: nucleus 



1. Introduction 

Near infrared observations of the central parsec of the Galaxy 
that were made over the past decade have brought new views 
of the environment in the vicinity of a supermassive black hole. 
They revealed a numerous population of young massive stars 
which may be distinguished into at least two different groups: 
within a distance < 0.03pc from the centre there are found more 
than ten so called S-stars orbiting the supermassive black hole 
on apparently randomly oriented and highly eccentric (<? > 0.8) 
orbits. They appear to be standard OB main sequence stars (Ghez 
et al. 2003, Eisenhauer et al. 2005) which is in contradiction 
with strong tidal forces that prevent stellar formation at the place. 
Unless these stars mimic their age, it is likely that they have 
migrated to the centre from larger distances. 

Further away, at 0.03pc < r < 0.5pc, nearly up to one hun- 
dred young stars were detected (see Paumard et al. 2006 for one 
of the most recent reviews). These stars are classified mainly 
as post main sequence OB supergiants and Wolf-Rayet stars. 
According to the evolutionary phase, their age is estimated to 
be 6 + 2Myr. Levin & Beloborodov (2003) pointed out that a 
substantial fraction of these stars form a coherently rotating disc 
(usually referred to as a 'clockwise' stellar disc or CWS). It is a 
flaring disc with an opening angle « 15° with a rather sharp in- 
ner edge at 0.03pc and it extends up to radius of ss 0.3pc. The ra- 
dial column density profile of the CWS decreases approximately 
as r~ 2 , i.e. most of the stars are concentrated at the inner edge. 
The mean plane of the disc can be determined by two angles: 
inclination i' x 127° with respect to the plane of the sky and 
longitude of the ascending node Q' ss 99° (measured from the 
North direction; see Paumard et al. 2006 for a detailed descrip- 
tion of the convention). Levin & Beloborodov (2003) suggested 
that this disc-like pattern is a consequence of a stellar formation 
in a self-gravitating accretion disc. 
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Further analyses (Genzel et al. 2003, Paumard et al. 2006) 
indicate the presence of another coherent stellar system which 
is usually referred to as the 'counter-clockwise' stellar disc 
(CCWS). This structure is assumed to be formed by fewer (< 
15) stars, it is narrower in the radial extent being concentrated 
around r > 0. 15pc and has a larger opening angle * 20°. The ex- 
istence of the CCWS disc is a matter of an ongoing debate (e.g. 
Lu et al. 2007), nevertheless, even if it is accepted as an explana- 
tion of the origin of another subset of young stars in the Galactic 
Centre, there would still remain more than twenty stars not be- 
longing to any of the two stellar discs and, therefore, without a 
satisfactory explanation of their origin. 

The gravitational potential in the considered region is domi- 
nated by the supermassive black hole of mass M. * 3.5xl0 6 M o 
(Ghez et al. 2003). It is surrounded by a roughly spherical clus- 
ter of late-type stars. The radial density profile is well fitted with 
a broken power-law p(r) cc with index (3 = 1.19 below 
r = 0.22pc and /3 = 1.75 above the break radius (Schodel et 
al. 2007). Its mass, M c within lpc is comparable to the mass of 
the black hole. 

The central region is surrounded by a molecular torus 
(circum-nuclear disc; CND) which lies at the outer edge of the 
black hole's sphere of influence (/?cnd ~ 1.5pc). Its mass esti- 
mated from the radio observations of ionised molecular gas is 
Mcnd ~ 1O 6 M (Christopher et al. 2005). This massive struc- 
ture defines a non-spherical component of the gravitational field 
in the central parsec. 

In this paper we investigate the influence of the CND upon 
the dynamical evolution of the disc-like stellar structures. In the 
subsequent section we briefly review the dynamics in the per- 
turbed Keplerian potential. In Sec. [3] we apply the results on the 
motion of the stellar discs — we present constraints on some 
parameters of the CND and the CWS determined from then- 
gravitational interaction and we also give suggestions on the 
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dynamics of the whole system of young stars over its lifetime. 
Conclusions and discussion of our results are given in Sec. [4] 



2. Dynamics in the perturbed Keplerian potential 

For the purpose of the study presented in this paper we intro- 
duce a simple model of the Galactic Centre which consists of 
three main constituents determining the gravitational field: (i) 
the central supermassive black hole of mass M. = 3.5 x 1O 6 M 
which is treated as a source of the Keplerian potential, (ii) the 
massive molecular torus modelled as an infinitesimally thin ring 
of radius /?cnd and mass Mcnd and (iii) a spherical stellar cusp 
with a power-law density profile (3 and mass M c within the ra- 
dius Rcnd- Both the ring and the cusp are centred on the black 
hole. 

The stars are treated as test particles whose motion is deter- 
mined by the composed smooth external potential. Their orbits 
can be represented by five orbital elements: the semi-major axis, 
a, eccentricity, <?, inclination, i, argument of pericentre, to, and 
longitude of the ascending node, Q. Here we assume the angles 
to be measured in the frame in which the ring lies in the x-y 
plane. For convenience, the results presented in the subsequent 
section will be transformed into coordinates with x'-y' represent- 
ing the plane of the sky and z' pointing from the observer. The 
inclination and longitude of the ascending node in this observer's 
coordinate system will be denoted i' and Q', respectively. 

If the gravity of the spherical cusp were ignored, the dynam- 
ics in the field of the central body and the ring would be equiva- 
lent to the reduced hierarchical three body problem. In this case, 
the orbital elements e, i, a> and Q will undergo secular evolution 
(Kozai 1962, Lidov 1962) on a time-scale of 



M. 



d3 
CND 



Mcnd a y/GM.a 
The equations of motion for mean orbital elements read: 

T K Vl - e 2 — =— e(l -e 2 ) sin2w sin 2 /, 
df 8 

Tu v 1 - e l — = e sin 2a> sin ; cos ; , 

df 8 

T K Vl -e 2 ^ = ^ {2 - 2e 2 + 5 sin 2 co [e 2 - sin' ; 

T K Vl -e 2 — =--cos/h + 4e 2 - 5e 2 cos 2 cj] 
df 4 L J 



]) 



(1) 

(2) 
(3) 
(4) 
(5) 



The temporal evolution does not depend upon the angle £2 which 
is merely a consequence of the axial symmetry; furthermore, en- 
ergy conservation implies a constant a in this order of approxi- 
mation. 

Including the gravity of the spherical cusp leads to an ad- 
ditional shift of the pericentre which can be incorporated by an 
extra term in equation @, dependent upon the global parameters 
of the cusp and the semi-major axis and eccentricity of the or- 
bit (Ivanov et al. 2005). The overall influence of the cusp can be 
characterised by a decrease of the amplitude of the oscillations of 
eccentricity and inclination and shortening of their period (Karas 
& Subr 2007). This is clearly seen also in Fig.Q]which shows the 
evolution of an example orbit both in the case with and without 
the potential of the spherical cusp. On the other hand, a generic 
influence of the cusp upon the evolution of Q. lies in diminishing 
the variations of its first time derivative; the characteristic time- 
scale of the change of Q. becomes in general much longer than 
that of the mutually coupled elements e, i and a). 
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Fig. 1. Evolution of orbital elements of two example orbits. The 
solid line represents a trajectory in the gravitational field of the 
central mass M. = 3.5 x 1O 6 M and a ring of radius Rqnd - 
1.5pc and mass Mcnd = M.. The dotted line shows an orbit 
integrated in the field including in addition a spherical cusp of 
mass M c = 0.1M.. In both cases the initial values of the orbital 
elements are: a = 0.1/?cnd, e = 0.1, i = 80° and Q. = 0. 



Within the context of this paper we are interested in a system 
where M c > O.IMcnd- In this case the amplitude of the oscilla- 
tions of eccentricity and inclination can be considered negligible 
and ai rotates with frequency much higher than that of Q. This 
configuration allows us to simplify eq. (0 by averaging over one 
revolution of a>: 



dQ 
~d7 



3 cos i 1 + 



3„2 



4 T K ^JT^e 2 



const . 



(6) 



The change of Q. over an interval Af can then be written as: 



AD. = — cos 1 a 
4 



= -5417° 
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According to the underlying perturbation theory, eqs. (Q]) - (0 
refer to the elements averaged over one orbital period, which 
cannot be trivially mapped to the osculating elements defined by 
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instant positions and velocities. We have performed numerical 
tests in order to check the error introduced by replacing averaged 
elements by the osculating ones in formula For M c ~ M, and 
e < 0.5 the difference between the analytical estimate and the 
numerically measured rate of change of Q was always smaller 
than a factor of 1.2. With increasing eccentricities, formula (O 
overestimates the real precession rate by a somewhat larger fac- 
tor which, however, still stays smaller than 2. 

The real source of the perturbing potential, CND, is rather 
clumpy torus of a finite thickness. Hence, we have performed 
several numerical integrations of a test particle orbit in a gravi- 
tational field of a central mass, spherical cusp and a set of ~ 10 
discrete point masses with orbits confined in a torus of a toroidal 
and poloidal radii 1 .5pc and 0.2pc, respectively. The difference 
in the orbital evolution with respect to that in the case of a 
ring-like perturbation was found to be negligible for orbits with 
a < 0.5pc. 

3. Consequences of differential precession 

3.1. Constraints on the CND 

Let us now consider an ensemble of stars forming a thin disc, 
i.e. with inclinations and longitudes of ascending nodes lying in 
a narrow interval. Let us further assume that the stellar orbits 
evolve solely due to the external gravitational potential deter- 
mined by fixed parameters M,, M c w M., Mcnd S M. and 
Rcnd S> The key feature of the orbital evolution will be pre- 
cession around the symmetry axis of the ring-like component of 
the gravitational field. 

If the semi-major axes of stars at the inner edge of the disc 
are smaller than those at the outer edge by a factor of > 5, the two 
edges of the disc will precess at a rate different by a factor > 10. 
After a certain period of time, their angular momenta will point 
to completely different directions, i.e. the disc-like structure will 
be destroyed. Hence, the requirement of the stability of the disc 
over a given period of time transforms into the requirement of a 
sufficiently slow precession at its outer edge. 

Let us consider the subset of the young stars in the Galactic 
Centre which form the 'clockwise' stellar disc (CWS). Inserting 
values Rcnd = l-5pc, a = 0.1/?cnd, Af = 6Myr and e — into 
eq. (0, we obtain 



AO. = -560° cos i 



Mcnd 
M. 



(8) 



In order to be compatible with observations, AQ has to be 
smaller than ~ 10°, which is the opening angle of the inner 
part of the CWS (Beloborodov et al. 2006). Hence, eq. ® 
poses constraint upon the inclination of the disc with respect 
to the molecular torus, depending on its mass. Considering e.g. 
Mcnd ~ 0.3M. (Christopher et al. 2005) requires cos i < 0.06, 
i.e. i e (86°, 90°). Simultaneously, this poses an upper limit 
6q < 5° on the initial opening angle of the stellar disc. For the 
sake of simplicity, we have considered a common sense of pre- 
cession of all disc stars, i.e. i < 90°. Identical results would be 
obtained for i > 90° due to the symmetry of the problem. The 
constraint would be tighter by a factor of a; 1.5 if we consider 
nonzero (~ 0.5) eccentricities of the stellar orbits. 

To conclude this analysis, we remark that cos i = news ■ 
iicnd < 0-06 is in accord with estimates of the normal vector^] 



1 We follow the convention of Paumard et al. (2006) according to 
which angles i' and SI' are related to the normal vector of the orbital 
plane as: (n' x , n',n') = (sin (' cos Si', - sin i' sin SI', - cos /'). 



of the plane of the disc and torus (n^ ws = (-0.12, -0.79, 0.60), 
Paumard et al. 2006; n' CND = (0.85, -0.40, -0.34), Jackson et 
al. 1993). 



3.2. A common origin of young stars in the Galactic Centre? 

Formula © indicates that orbits of stars at radii > O.lpc and/or 
inclinations i < 85° or i > 95° were considerably affected by 
precession within the past 6Myr, i.e their current orbital param- 
eters are different from their values at the time of the birth. We 
suggest a possibility that stars which are not considered to be 
members of the CWS nowadays have been its members at the 
time of its formation. During the ~ 6Myr of the dynamical evo- 
lution their orbits were subject to precession due to the gravity 
of the CND and were detached from their parent stellar system. 
This model could represent a possible solution of the problem of 
the origin of all young stars in the Galactic Centre. 

The mapping between the initial and the current orientation 
of the stellar orbit is formally rather straightforward within our 
simple model. Unfortunately, the observational data do not pro- 
vide us with sufficiently accurate values of the parameters Rcnd 
and Mcnd- Furthermore, the high sensitivity of the precession 
rate upon the inclination and semi-major axis together with a 
lack of robust determination of these orbital elements from the 
observational data also stand as a severe obstacle in an attempt to 
track the orbits of the observed stars back in time, which could 
prove or discard the hypothesis of a common origin. In the rest 
of this section we describe a test shows that our model is com- 
patible with the publicly available observational data. 

We have taken data from Table 2 of Paumard et al. (2006) 
from which we have considered all stars with determined 3D ve- 
locity and index > 15 (i.e. excluding the S-stars) which gives in 
total Af* = 72 stars. Five free parameters of the model consist of 
the two angles, (Q.' , z'l), determining the initial orientation, nl, 
of the stellar disc; another two angles, (Q CND , i' CND ), determine 
the orientation of the CND, and Mcnd represents its mass. (The 
last parameter can be considered as a degenerate combination of 
Mcnd, ^cnd and Af; in the following, we will implicitly assume 
Rcnd = 1 -5pc and At = 6Myr.) For a given set of parameters we 
scan the 1 cr neighbourhood of each star's velocity with sampling 
d v . The x? and y' coordinates of the stars' positions are assumed 
to be determined exactly. On the other hand, the z' coordinate 
(along the line of sight) is unknown. Therefore, we scan it with 
sampling d z in a whole range allowed by the condition that the 
star is gravitationally bound to the black hole. In total, we con- 
sider Vicr = d z dl pairs of position and velocity vectors which 
represent states compatible with the observational data of a par- 
ticular star. For each state we perform a rotation of the normal 
vector of the orbit around the axis of the CND according to for- 
mula (O, which gives its direction, n^o at f = 0, i.e. 6Myr ago. 
We further calculate its angular distance to no, cos 6q = n 7i o ■ no 
and count the number of states, Nj$, with So < 5°. We con- 
sider the measured star's position and velocity to be compatible 
with the hypothesis that it was born in the disc with normal vec- 
tor n' Q and thickness 5°, provided Np > 0. Finally, we denote 
Ns(Q' , i' , i^cND' 'cnd' -^cnd) the number of stars with Nj t s > 
for a given set of values of the parameters of the model. 

As we have discussed in Section 13.11 the inner part of the 
CWS must have undergone only negligible precession due to 
the gravity of the CND. Therefore, we assume that it conserves 
orientation of the putative single parent disc and we consider 
(Q.' Q , i' ) = (99°, 120°) which is the normal vector of the inner 
part of the CWS according to Beloborodov et al. (2006). We 
further set Mcnd = M. which enables us to plot N5 as a func- 
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Fig. 2. Various tests of compatibility of the hypothesis of a common origin of young stars in a single thin disc. N5: number of stars 
that have at least one state in the Icr neighbourhood of observed velocities that corresponds to an orbit dragged from the parent disc 
with normal vector (Q ( ' ( , i' Q ) — (99°, 120°) by the CND of mass Mcnd = M.. N e ^: similar to the previous but excluding orbits with 
e > 0.5. Vs. volume of the subset of the lcr neighbourhood occupied by the states originating in the parent disc. (Ms): mean value 
of the mean anomaly of orbits with 60 < 5°. 



tion of Q' CND and /p ND as it is shown in Figure [2] We see that 
there exists an extended region where the observational data of 
nearly all stars are compatible with the hypothesis of their origin 
in a parent disc of thickness ~ 5°. This region of large values 
of N5 extends along the set of (£2 CND , i' CND ) perpendicular to the 
normal vector n ( ' ( . This is a natural consequence of the assump- 
tion that ~ 35 stars, identified as CWS nowadays, haven't un- 
dergone large precession. The region of good compatibility also 
includes an approximate orientation of the CND as determined 
from observations, (Q^ ND , i' CND ) = (25°, 70°), e.g. by Jackson et 
al. (1993). 

We have performed analogical test of compatibility also for 
Mcnd = 0.3M. and 3M. and ~ 10° neighbourhood of (Q.[ r i' Q ) = 
(99°, 120°). In all cases we have obtained a picture similar in 
that there exists an extended region with N5 > 68. Enlarging the 
inspected neighbourhood of the observed velocities to 3cr leads 
to larger values of N5 with its maximum reaching 72. This means 
that the observational data are compatible with the hypothesis of 
the common origin of the young stars in a single thin disc, never- 
theless, they do not pose strong constraints on the parameters of 
the model. We introduce three supplementary tests that may be 
applied to the observational data to verify validity of our hypoth- 
esis. First, the model of star formation is assumed to prefer low 
eccentricities of the stellar orbits. Hence, we introduce N e ^ in the 
same way as N5, but now with an additional condition e < 0.5 
for all tested orbits. Middle panel of Fig. [2] shows that region 



with large N e> s coincides with the region of large N5. Maximum 
value of N e ,5 is ~ 55, i.e the model requires about one third of 
the stellar orbits to have moderate to large eccentricities. Again, 
considering 3cr neighbourhood of the measured velocity vectors 
weakens this constraint, giving N e .s > 70 for a wide range of the 
model parameters. Additional analysis reveals that most of the 
eccentric orbits does not belong to the CWS subset of stars. 

A ratio Nj^/Via- can be considered as a measure of the prob- 
ability that the orbit of star j originated in a disc with opening 
angle 5° and normal vector w n . Consequently, we introduce 



r 5 = 



(»* N + 1^* 



(9) 



as a measure of the probability that all stars originated in a 
thin disc (Adding a unity to Nj t $ in (0 prevents P5 from be- 
ing zero everywhere while it does not strongly affect its mean- 
ing.) Bottom panel of Fig. [2] shows that ^5 accentuates orien- 
tations of the CND nearly perpendicular to the normal vector 
(Q.[ v i' ), but it does not strongly discriminate among the mod- 
els that fall into this region. Maximum of this function is at 

("cnd> <cnd) * (55°, 40°). 

Following Beloborodov et al. (2006), the configurations 
compatible with the hypothesis of a single warped disc are ex- 
pected to have equally distributed value of the mean anomaly, 
M, of the individual orbits. Full test of the distribution of the 
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mean anomaly of all configurations that have fulfilled other cri- 
teria of compatibility is not possible as it would require analysis 
of « PsV^* combinations of orbits. (Note, that other tests pre- 
sented here require analysis of only N*V\, a individual orbits.) 
In the bottom panel of Fig. |2]we present a restricted test show- 
ing the mean value of the mean anomaly, (M5), for all tested 
orbits with d>o < 5°. This quantity is close to n which corre- 
sponds to the expected uniform distribution in the major part of 
the (fi^ ND , i' CND ) space. Analogical plot would show that (M%) 
is close to the expected value ^n 2 in the regions of large value 
of /V 5 and V 5 . This indicates that our hypothesis does not require 
some preferred value of M and, therefore, some of the configu- 
rations with large N5 are also compatible with the assumption of 
random distribution of the orbital phases. 



4. Conclusions 

The massive molecular torus (CND) surrounding the central par- 
sec of the Galactic Centre causes precession of the orbits of 
young stars which move at distances 0.03pc < r < 0.3pc around 
the supermassive black hole. The rate of the precession depends 
on the orbital parameters as well as on the orientation and mass 
of the CND. This rate is comparable to the lifetime of the young 
stars for a wide range of parameters and, therefore, this process 
should be taken into consideration in attempts to determine the 
relation between initial and current values of their orbital pa- 
rameters. We have shown that Mcnd Z 0.3M. would destroy 
any coherently rotating disc-like stellar structure within 6Myr, 
provided the inclination of most of the orbits with respect to the 
CND deviates by more than 5° from 90°. In other words, the sta- 
bility of the stellar disc within its lifetime poses constraints on 
its inclination with respect to the CND and on the mass of the 
CND. 

We further suggest that most if not all young stars observed 
in the Galactic Centre may have been formed in a single, initially 
coherently rotating structure, presumably via fragmentation of 
a thin self-gravitating gaseous disc. Within this hypothesis, the 
orientation of the stellar disc was nearly exactly perpendicular 
with respect to the CND. Its 'core', represented by the CWS 
nowadays, remained nearly untouched by the precession. On the 
other hand, stars that were formed at the outer parts of the disc 
and/or slightly off its mean plane, or that were scattered out of it 
via two-body encounters, have undergone a more rapid preces- 
sion of their orbits, i.e. they apparently don't belong to the stellar 
disc any longer. We have shown that within the 1 cr uncertainty of 
their current velocities there exist such parameters of the stellar 
orbits that would have had their angular momenta collinear about 
6Myr ago. Due to the high sensitivity of the precession upon the 
orbit inclination with respect to the CND and the uncertainty in 
the observed parameters of the stellar orbits, the procedure de- 
scribed in the previous Section cannot provide robust constraints 
on the parameters of the model. Therefore, the concept of a sin- 
gle warped disc of young stars in the Galactic Centre may be 
considered as being viable, but not proven yet. Our hypothe- 
sis, however, gives an explicit prediction on a specific pattern 
of the normal vectors of the stellar orbits which may be deter- 
mined from future, more accurate observations: all of them are 
assumed to be found close to the circumference perpendicular to 
the normal vector of the CND. 

Let us emphasise that the gravitational influence of the CND 
leaves stronger imprints on the dynamics of stars more distant 
from the centre. Hence, we suggest that these stars deserve fur- 
ther attention from the observational point of view. Improved 



measurements of their kinematical state may bring a new light 
into the question of the formation and dynamical evolution of 
the population of the young stars in the Galactic Centre. Beside a 
generic demand on better constraints on orbital parameters from 
the observational side, there is also a room for improvements of 
the model itself. Its most important (and computationally rather 
expensive) modification will probably lie in an improved treat- 
ment of the evolution of the individual orbits, which would take 
into account gravity of the stellar disc itself. 

As a final remark let us note that the strict constraints on 
the mutual (perpendicular) orientation of the stellar disc and the 
CND raises a question about the dynamics of gas from which 
the young stars were formed. It is likely that the parent gaseous 
disc had to be nearly perpendicular to the CND, so that it would 
not be destroyed via differential precession before it gave birth 
to the numerous stellar population. Such an initial orientation 
is statistically not very probable, opening the question whether 
it can be a generic result of dissipative (hydro)dynamics in the 
resonant external potential. 
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